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Abstract 


Theoretical research has been carried out to study the effect of free-stream 
turbulence on sonic boom pressure fields. A new transonic small-disturbance 
model to analyze the interactions of random disturbances with a weak shock 
has been developed. The model equation has an extended form of the clas- 
sic small-disturbance equation for unsteady transonic aerodynamics. An 
alternative approach shows that the pressure field may be described by an 
equation that has an extended form of the classic nonlinear acoustics equa- 
tion that describes the propagation of sound beams with narrow angular 
spectrum. The model shows that diffraction effects, nonlinear steepening 
effects, focusing and caustic effects and random induced vorticity fluctu- 
ations interact simultaneously to determine the development of the shock 
wave in space and time and the pressure field behind it. A finite-difference 
algorithm to solve the mixed- type elliptic-hyperbolic flows around the shock 
wave has also been developed. Numerical calculations of shock wave interac- 
tions with various deterministic and random fluctuations will be presented 
in a future report. 
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1 Introduction 


Experimental data exists showing the pressure profiles of sonic booms, cre- 
ated by the passage of a distant supersonic aircraft, can be drastically 
affected by free-stream atmospheric turbulence. 1-5 The usual /V-wave or 
a shaped sonic boom profile can randomly exhibit either large pressure 
peaks with short rise times or rounded profiles with longer rise times or 
messy pressure signatures. Recent laboratory model experiments to study 
the effect of turbulence on the rise time and wave form of N - waves have 
shown similar results. 6 The interaction of the sonic boom with the atmo- 
spheric turbulence, specifically in the atmospheric boundary layer relatively 
near the ground, may result sometimes in higher, and may be unacceptable 
loudness levels.' Therefore, in order to get reasonable estimates of the sonic 
boom performance of various designs of a future supersonic transport air- 
plane it is essential to understand the basic interactions of the atmospheric 
turbulence with shock waves. 

The basic analysis of the distortion of sonic bangs by atmospheric tur- 
bulence was given by Crow. 8 Using a first-order acoustic scattering theory, 
Crow showed that the pressure perturbation behind the shock is related to 
the interaction of the shock with the disturbances it encounters while mov- 
ing in the atmosphere. The pressure profile can be calculated by a surface 
integral over a paraboloid of dependence, whose focus is the observation 
point and whose directix is the shock front. By describing the turbulent 
eddies in the Kolmogorov inertial subrange, it was found that the mean- 
square pressure perturbation behind the shock changes like (A p) 2 (t c /t) 7 ^ 6 
where (A p) is the pressure jump across the shock, (t) is time after the shock 
passes an observation point and (t c ) is a critical time predicted in terms 
of meteorological conditions. Crow’s analysis predicts reasonable average 
values of the pressure fluctuations for times ( t ) comparable to (t c ). 

The singularity in the pressure perturbations near the shock front (when 
( -* 0) was analyzed by Plotkin and George 9 . A second-order acoustic 
scattering theory was used to describe shock rounded signatures. The av- 
erage of the diffraction effects was approximated as a dissipation term. 
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The thickening of the shock is explained as a balance between nonlinear 
steepening effects and the dissipative effect of the turbulent scattering of 
acoustic energy out of the incident shock. Rise time predictions of this 
theory show some correlation with experimental data. On the other hand, 
Ffowcs Williams and Howe 10 examined the approaches that describe the 
possibility of a turbulent thickening of weak shock waves and reached a 
conclusion that atmospheric turbulence cannot be the cause of shock thick- 
ening. They suggested that weak shocks may attain a dispersed profile due 
to non-equilibrium gas effects. 

It should be emphasized here that the scattering analyses of References 
8-10 considered small turbulent perturbations against the shock strength, 
whereas in the case of the interaction of the sonic boom with atmospheric 
turbulence the flow random fluctuations may be of the same order of the 
shock weak strength and may strongly distort the shock front. The above 
analyses also did not account for shock jump conditions that must be sat- 
isfied in an inviscid analysis across any shock surface. The approximation 
made in Reference 9 of the average diffraction effects described as a dissi- 
pative term is also unclear. 

A different approach was taken by Pierce 1112 . He interpreted the spikes 
observed on sonic boom pressure waveforms as being due to the simultane- 
ous focusing and diffraction of a nearly planar .V-wave by an inhomogeneous 
layer in the atmosphere. The shock front develops ripples that are trans- 
formed into folds in the front when the shock passes vertices of caustics. 
This mechanism results in a fine structure of very small pressure jumps 
that correspond to the various segments of the folded wavefront. Pierce 12 
derived a stochastic model of a sharp shock propagating through a turbu- 
lent atmosphere to substantiate the very small discrete structure of sonic 
boom profiles. 

Pierce model 11 ' 12 , however, neglects nonlinear effects that become sig- 
nificant specifically near a caustic vertex as was shown by Cramer and 
Seebass 13 and Gill and Seebass 14 . Cramer and Seebass 13 described the fo- 
cusing of a very weak and slightly concave shock wave by the unsteady 
transonic small-disturbance flow equation. Gill and Seebass 14 derived an 
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approximate analytical solution of the steady transonic small-disturbance 
problem for the nonlinear behavior of a weak compression wave with a finite 
rise time near a caustic. They calculated the reflected shock wave from a 
caustic and provided an estimate of its strength. The experimental results 
of Sturtevant and Kulkarny 15 also show that focusing effects are specifically 
significant for weak shock waves as occurs in the case of the sonic boom 
signatures. 

Sparrow and Pierce 15 have recently presented a simple statistical predic- 
tion for how often sonic booms propagating in the earth’s boundary layer 
will encounter caustics. The theory is based on describing the variation 
of ray tube areas of a sound wave propagating in a turbulent medium by 
a generic harmonic oscillator equation. For realistic realizations of atmo- 
spheric turbulence the model predicts that sonic booms will exhibit spikes 
with the occurrence of caustics after a very short distance of travel in the 
random medium, thus agreeing with the predictions of Pierce 11,12 . 

In a recent paper. Pierce 17 has derived a model equation to describe 
the development of sonic boom signatures in an atmospheric turbulence. 
The equation has been constructed by suing logical physical considerations 
only. It extends geometrical acoustic approximations to include convec- 
tion at the wave speed, diffraction effects, molecular relaxation, classical 
dissipation and nonlinear steepening effects. The atmospheric turbulence 
enters through an effective speed of sound which varies randomly in time 
and space. However, since this theory has not been developed consistently 
from the fluid dynamic equations, Pierce 17 raised questions whether ail the 
effects are necessary in his suggested model and how to accomplish a nu- 
merical or analytical solution to the problem. 

Related with the problem of the sonic boom interaction with turbulence 
is the basic question of the interaction of a shock wave with a vortical flow 
and specifically with a vortex or a train of eddies 18-21 . The research of 
the later problem was basically motivated by the interest to reduce the 
noise and vibrations produced by high-speed supersonic vehicles. For these 
problems, the interaction of relatively strong shock waves with turbulent 
jets or wakes is a significant source of noise. The shock-vortex system is a 
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basic element of these more complex interactions 21 . It can also shed light on 
the sonic boom interaction with atmospheric turbulence, specifically when 
the shocks are weak anu vortex strength is comparable with the shock 
strength. 

Experimental results of shock-vortex systems 1 ® -21 revealed curved, dif- 
fracted shocks as well as complicate structures of reflected curved shock 
waves from the incident shock front due to the vortex induced flow field. 
The pressure field behind the shock contains regions of compressions and 
rarefactions that produce acoustic waves. Similar shock structures were 
also observed by Sturtevant and Kulkarny 15 who investigated the focusing 
of weak curved shock waves. Of specific interest are Dosanjh and Weeks 19 
measurements of the interaction of a shock wave with a vortex street. The 
shock front is distorted by the wake flow which probably results in a focusing 
process, while the vortex street is rapidly dissolved by the shock. 

The analyses of the shock-vortex interaction are limited to linear per- 
turbation theories only . 22 26 These analyses considered the jump conditions 
across a shock surface and predicted the development of vorticity waves, 
entropy waves and acoustic waves behind the shock front . 22-2 * 1 The acous- 
tic wave was approximated by a quadrupolar 26 or as a sum of monopole, 
dipole and quadrupole acoustic sources. However, since sill of these theo- 
ries are linear, they cannot account for any nonlinear effects due to shock 
large distortions, focusing and caustic effects or nonlinear steepening ef- 
fects that are found in experiments 18-20 or in recent numerical simulations 
of shock-vortex interactions . 27-29 

The review of experimental and theoretical investigations of the interac- 
tion of shock waves with free stream vortical or turbulent flows shows that 
this complex nonlinear interaction is still an open problem . 27 Specifically, 
the improved simulation of sonic boom propagation through the real at- 
mosphere requires a better understanding of the interaction of sonic boom 
with atmospheric turbulence . 16 

The analysis of the experimental data and the theoretical approaches 
shows that in the case of the sonic boom, the shock waves near ground 
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are very weak, but still stronger than any acoustic wave. Also, flow fluc- 
tuations due the atmospheric turbulence can become comparable to the 
shock weak strength such that locally the shock strength can be strongly 
reduced or magnified and the shock wave front can be distorted significantly. 
Therefore, linearized acoustics and its second-order scattering problem, or 
first-order linear theories of shock- vorticity interaction do not represent cor- 
rectly the development of the weak shock and the pressure field behind it 
(see also Section 2). However, in a coordinate system moving with the basic 
weak shock, the problem may fit the transonic framework. 

This paper presents a new transonic small-disturbance model that has 
been developed to describe the interactions of random fluctuations with a 
weak shock wave. The model equation is found to have an extended form 
of the classic nonlinear acoustics equation that describes the propagation 
of sound beams with narrow angular spectrum (KKZ equation). 30,31 The 
model shows that diffraction effects, nonlinear steepening effects, focus- 
ing and caustic effects and random induced vorticity fluctuations interact 
simultaneously to determine the development of the shock wave in space 
and time and the pressure field behind it. A finite- difference algorithm to 
solve the mixed-type elliptic-hyperbolic flows around the shock wave is also 
presented. The results of the numerical calculations will be presented in a 
future report. It is expected to find solutions that will describe both peaked 
or rounded or messy pressure signatures as were recorded in experiments. 


2 The Breakdown of the Linearized Theory 

2.1 The Linearized Theory 

An inviscid and a non heat conducting flow is assumed. A normad shock 
with a uniform supersonic oncoming stream and a uniform subsonic outgo- 
ing flow is considered. The upstream flow ahead of he shock is characterized 
by a speed ( Uo a ) in the x- direction, pressure (p 0a ) and density (p 0a ) and 
the downstream flow behind the shock by (t' 0 <>), (pob) and (pob) respectively. 
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Assuming the shock front is given by the x = 0 plane, the jump conditions 
across a normal shock 32 show that: 

POat'Qa = PobUob 

POa^Oa+POa = Pob^Qb + Pob (1) 


1 „ 7-J-3 , 7 _ rr 

POa t'oo ■" POa Oa 

2 7-1 

Small disturbances are considered 


_ 1 7t 3 , 7 rr 

— ~^POb '- Ob -po6C'Of> 

- 7-1 

in each of the uniform streams. 


The 


velocity vector (V). pressure ( P), density {p j and vorticity (w) are given 
ahead (j = a) and behind the shock (j = 6) by: 


V = i'o]{e +cu H ) 

- } ~r 

Z 3 , = Pq } { 1 + + • ■ ■) 


( 2 ) 


Pj = Poj(l + tp\ } + • --) 


= ew + • ■ • , w .= y x u 

~ 1; ~ lj -1; 


Here v ,p 1; , p 1; are functions of (x,y,z,t). An axial coordinate moving 

with the uniform speed is considered in each region, — x — U 0) t, T = 
U 0 jt. The substitution into the continuity, momentum and energy equations 
results to the leading order in (V = -h-r-e +^-e): 

< - V~ y ■>- -j 
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Equations (3) result in: 


<W 

dr 


+ 7(V-i’ ) = 


d P\] dpi } 

a7 = 7 ^T 


du 

-Uj. 

dr 


= 0. 



d 2 Pi j 
dr 2 


- V 2 Pu = 0 


(4) 


where M 0} = T 0; /c 0j and c\ } = 7 Poj/poj- Equations (4) show that u> = 

2 ) an< ^ t ^ lat t ^ ie pressure perturbation p X] is described by the acous- 
tics equation. Therefore, the first order disturbance flow can be split into 
a linear combination of rotational and irrotational parts: v = v + v 

The rotational part can be described essentially by incompressible flow 
equations: 

Y V ~^ = 0 ' (5a) 

The irrotational (potential) part may be described by acoustics equations 
relative to the basic flow in each region; 


v 

* ijo 



_ , r2 d&j 

Pl) ~ 


M, 


2 d 2 <t>, 


Oj 


dr 2 


- V 2 <£; = 0 


(5b) 


The first-order perturbation theory also considers the distortion of the 
shock front. Assuming that the perturbed shock front is given by x — 
= 0 , the exact jump condition across the shock 32 result to the 
leading order in a set of conditions that must be satisfied along the x = 0 
plane for any (y,z,t): 

MO, y, *, t) + MO, y, s, f ) - £ 1 = P\b( 0. y. t) + u 16 (0, y, z, t) - ill 

U 0a U 0b 


U Qa ^Pio(0, y,z, t) + 2u la (0, y, s,t)) + -1—^(0, y, z, t)J = 

= Un ^pi6(0, y, r, t) + 2ui(,(0. y, t ) + — ^ p^O, y, r, t)^ 
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U()a Pla(0* y ' + * t) + 3u\ a (0) y. Z.t) + 


(7 - DXIL 


■(P\ a {0,y,s,t) + 


+ «i«(0,y,z,O) - ~ I 1 + 


2/7 


U 0a V (7 — l)iV/(ja, 


= £’ot ( /W 0 - y 1 s.f) + 3 u, 6 (0. y.r.f ) + 


(7 - l)-V/ 0 \ 


(Pi»(0, y, z,t)+ 


4- u lt (0. y. -.<)) - ^r 1 (l + T ^Ty\77 

lob\ (7-l)-Vf 0 J 6 


^Oo(t f io(0, y, r,i) + y ly ) = To6(i>i&( 0. y, r, f ) + y ly ) 


koa( ^ia(0, y, 0 + g\z) = L'o 6( u’ii(0, y, ~, 0 4- yi, ) (6) 


Here y i( = dgi/dt and (ui ; , u'i, ) axe the components of the velocity 
perturbation . The linearized jump conditions in Eqs. (6) include the 


entropy increase produced by the shock. It can be shown from Eqs. (5) and 
(6) that, using the solution of the downstream equations for the disturbance 
flow, the shock conditions are adequate to describe the flow downstream and 
the disturbed motion of the shock wave for given upstream disturbances. 


2.2 One Dimensional Flow 


In the case of a one-dimensional flow, the rotational part vanishes identically 
and the solution of the acoustics equation ahead (j = a) and behind ( j = b) 
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the shock is given by: 



P„ - Mo, ,J 0) ) - G 

where F and G axe arbitrary functions that describe the upstream and 
downstream acoustic waves. Here p 3 is an arbitrary function that describes 
entropy waves that axe convected with the flows (the entropy first-order 
disturbance is given by 5 ; - 5 0j = = —ec p M 0j p 3J (£ j), where S } is 

the entropy in region S 0j = c v ln(P 0j /Po; ) and c„,c p are the specific 
heat constants). Assuming that no upstream acoustic waves can develop, 
specifically not in the flow behind the shock then G = 0. Then the shock 
jump conditions (6) provide a system of 3 linear equations for the solution 
of the downstream acoustic and entropy waves F t (£(, - and p 3 t>{£b) and 
the shock position rate of change in time <7 k(*) in terms of the given acoustic 
and entropy perturbations F a (£ a - and p, a (<f a ) in the upstream flow. 
Let, 


F*0 F a ( ( (~ 0a Q)a ) 0 1 PsaQ — Psa{ f^OaO 


( 8 ) 


FbO Fb{ ( "1“ Cot)^)- PsbO “ P$b{ 


then from Eqs. (6) - (S) we get: 
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1 4 Mob)FbO 4 MobpsbO 4 f 77 l) J7~ = (1 4 -V/oa)-f’aO 4 M OaPsaO 

VtOa ' <- Ob 


^2 4 Mot, + F* 0 + MobPsbo = ^2 4 M 0a 4 'jFao 4 M 0a p, a0 


3 4- Mot, 4 - 2 ( 1 4 7 ^ 04 ) ) P'bO 4 Mobpsbo 4 

v (7 - 1)A/q6 / 


+(fe(l+ . 2 /- 

U r 06 V (7-l)-V/ 0 2 a 


, 2/7 \ gn _ 
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The determinant of the system (9) may be written in the form A = 
(M$ a — l)fn(Ml). Therefore, the solution of Eqs. (9) shows that the 
shock front motion and the perturbed flow behind it may be described by: 


X as 


ML - 1 


9 \( t\ Mq<x ) 
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ML - 1 


A/ Oo ) 4 
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Plb — Pob 1 4 
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(10) 


The functions ^i,u ll pi,p l can be expressed in terms of the given flow per- 
turbations ahead of the shock wave. In principle, these expressions may 
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enable a spectral characterization of the pressure fluctuations and the tur- 
bulence downstream of the shock wave in terms of the spectral character- 
ization of the incoming turbulence. However, Eqs. (10) show that the 
linear approach is a non-uniform approach when the shock wave is very 
weak M$ a — ► 1* and the flow fluctuations are of the same order of the 
shock strength, e ~ — 1), as is the case of the interaction of the sonic 

boom with atmospheric turbulence. A similar nonuniformity is also ex- 
pected from the analysis of two- or three-dimensional flows. However, this 
uniformity problem leads to a different approach to study the interaction 
of a weak shock with comparable random fluctuations in the flow. 


3 A Transonic Small Disturbance Model 


The analysis of the linearized problem of the interaction of a weak shock 
with small disturbances shows that it is an invalid approach when the flow 
perturbations are of the order of the shock strength. Therefore, a different 
approach has been developed to study the interaction of weak shocks with 
comparable random fluctuations in the flow. In a coordinate system moving 
with a basic given weak shock, the problem may fit the transonic theory 
framework. A transonic small-disturbance model is developed to analyze 
the flow across a basic weak shock running in the ( — x) direction. A coordi- 
nate system attached to the basic shock is considered. The velocity vector 
pressure (P), density ( p ) and vorticitv (-,•) are described every where 
in the flow by: 

V = C/oo{t'(l+€ 2 / 3 U+CU l +e 4 / 3 U 2 + - -) 

+;(«h + e* P v 2 H ) + k(eWi + e* / 3 w 2 + •■•)} 

P = Poo (l + e 2/3 P + e Pl + e 4 / 3 p 2 + •••) (11) 

P = P> o( 1 + e 2 ^ 3 P + e P i 4- e 4 ^ P j 
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j Wy + kw z ) + t{ luJzi + j^y\ + k'^'zi ) + • ■ • 

where U » = a*(l + fe 2/3 ) is the speed of the basic shock {I< >0) and 
^OO > Poo » Poo are the speed of sound, pressure and density of the unperturbed 
flow ahead of the shock, (e 2 ^ 3 ) represents the scale of strength of the basic 
weak shock where e 1. A rescaling of the x-coordinate and time ( t ) 
has also been considered: x* = and t m = fa x .« l/3 , such that each of 
the terms in (11) is a function of (x\y, s,f’). The rescaling in x means 
a stretching of the picture of the flow around the basic shock in order to 
capture the basic nonlinear effects that occur in the flow across the shock. 
The rescaling in time accounts for low-frequency unsteady perturbations in 
the flow. The constant A reflects that the speed of the basic shock wave is a 
little higher than the speed of sound, ahead of the shock. The substitution 
of Eqs. (11) into the continuity, momentum and energy equations results 
to the leading orders in: 


dx‘ 


( P 4" U ) as 0 
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dx 


;(Pi + u,) = 0 


dp d dv x 

W + a? {Ul + P2+pu) + -^ + 


dw\ 

~d7 


= 0 


~;(7U + p) = 0 

1 £:<vu+p,) = o 


du dui du 1 dp -2 

W- + gp + <' “ + 0 + >' 1 gp + ~gp - 0 


( 12 ) 


(13) 
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dv\ 
dx‘ 

du.'\ 

( dx m 


dp 

df 


dp dp 


1 dp 

7 dy 

1 dp 
7 dz 

dp dpi dp2 




d. r* <5x ' 


dx* 


( 14 ) 


(15) 


d W* 


= 0, 


dw- 


= 0 


(16) 


<9x* 5 x" 

From the equation of state and the definition of entropy it can be shown 
that the temperature T m and entropy 5 are given by: 


where 


r m = t^( \ + t 2/3 T + -■■) 

S = S.y^i 1 4 " + • 1 ■ ) 


T» = p^/Rp^. S~c = c„ In ( ^ 

\^°°/ 


C 

T = p - p, S = -~{p - 7/5) 


Equations (12) through (17) result in: 

u + P = /(</• r.f) 


(17a) 


(17b) 


(18a) 


7 U +p = y(y,c,t*) 


(18b) 




di't dw\ 


- 2— +{-K - f + g -( 7 + l)u)— + 7 J 7 + 


< 9 x* dy dr 


1 dg 

(,8c) 
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y 


(13d) 


du div i _ 1 dg 

dz dx * 7 5c 

du dt’i _ 1 dg 

dy dx m 7 dy 


(18e) 


r * - / - (7 - l)u, ■i = ~{g-yf) (18f) 

X 

where / and y are random induced fluctuations due to the free turbulence. 
The function g is related to the vorticity fluctuations in the flow. Equations 
(18) show that the axial perturbation (u), pressure perturbation (p) and 
density perturbation (p), that are of order of the shock strength (e 2 ^ 3 ), 
interact with the transverse velocity perturbations c, and uq. that are of a 
smaller scale (e). 


The substitution of u — g/~ ; + u in (lSc). (ISd) and (18e) results in a 
problem for solving a velocity potential function o(x’. y. z, t' ) where: 




do 

dx m 


~ do 

dy' 




(19a) 


-<t>i 


•t* 



- + /+( 7 + DOr* 
7 


^r*r* ( Oyy 4“ d • 2 ) 


In a conservative form Eq. (19b) is given by: 


1 dg 
7 dt m 


(19b) 


^2 4> x > + — + ((A 4- gj‘t + f)o I • + (7 + l)<j> 2 ./2)) x « — (o y ) y — (<t> z ) z = 0 

(19c) 

The exact shock jump conditions (Ref. 32) must be satisfied along any 
shock surface x m — h(y,z,t’) = 0 that may appear in the solution. To the 
leading orders they result in: 

(/] = 0 . [y] = 0 ( 20 a) 
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“ 2 [^*l^ + ( A ’ + 5 + /)[<fi*] + (7 + l) = 0 (20b) 

r i r i dh . . _ dh 

K1 + [<M^=° [<t>z] + {<t>Aj- z = 0- (20c) 

where [a] represents the jump across the shock property a, [a] = a B — a A . 
Equations (20a) show that to the leading order there is no jump in entropy 
across the shock, [5] = 0. Equations ( 11 ) and ( IS) also show that the local 
Mach number Me at any point in the flow is given by: 

M? - l = e 2/3 u m . u' = |(7 + l)Ox- +A' + /+ (21) 

The flow is locally supersonic when (7 -f l)<£ r . -f A' -+- / + * >0, sonic when 
(7 -f 1 )<Pz- + A’ 4- / + * = 0. and subsonic when (7 + 1 )<p x » + K + / + * < 0. 
Equations (19) and (20) are an extended version of the classic unsteady 
small-disturbance equation for transonic aerodynamics (Cole and Cook 33 ). 
The only changes are due to the random terms g and /. Starting from 
given functions for / and g and initial conditions that describe a given 
basic shock, Eqs. (19) and (20) can be integrated in space and time to 
describe the development of the shock wave and pressure field behind it. A 
numerical algorithm to solve these equations is described in Section 4. 

An alternative approach may be found by taking a x m - derivative of ( 18c ) 
and using Eqs. (18a) and (14). The pressure perturbation (p) satisfies the 
equation: 

d (dp K + / + g_h 7 + 1 dp_\ _ 1 ( &J> &Y\ 

dx- + 2 27 P dx -J 2 \dy 2 + ds 2 ) (2 “ } 

Equation (22) is an extended version of the classic KKZ equation that 
describes the propagation of nonlinear sound beams with narrow angular 
spectrum in an inviscid fluid (Zabolotskaya et al. 30 . Kuznetsov 31 ). Eq. 
(22) also has a similar form to the model equation that has been recently 
developed by Pierce 1 ' using logical considerations only. 
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Equations (19) and ( 22 ) show that diffraction effects, nonlinear steep- 
ening, focusing and caustic effects and random induced fluctuations due 
to turbulence interact simultaneously to determine the development of the 
shock wave in space and time and the pressure field behind it. Turbulence 
tends to change the local speed of sound in the flow across the shock and 
through this effect to reduce or to magnify the strength of the jump along 
the basic shock (see Eq. (21)) or to distort the shock front. These changes 
may result in unsteady motion of the shock front or in caustic vertices or 
in reflected shocks behind the incident wave that can produce the variety 
of pressure signatures of sonic booms that are measured in experiments. 


4 Finite Difference Scheme 


A finite difference algorithm to solve the unsteady mixed- type elliptic- 
hyperbolic flow around the shock wave has been developed. Murman and 
Cole 34 and Cole and Cook 33 techniques are used. A fully conservative 
scheme that is based on the conservative form (Eq. 19c) is used. In this 
way the difference equations also contain the shock relations (Eqs. ( 20 )). 


Consider a uniform finite difference mesh (Ax*. Ay, Az, At*) in space 
and time, with points (x*,y,c,f) labeled by ( i.j.k,n ). The results can be 
easily generalized to a variable mesh. Eq. (19c) can be expressed in a con- 
servative flux form for a box centered on a mesh point ( i,j,k ). Therefore, 




A + —g 4- f J <p r - + (7 + 1 )<p] 


' (»+£,;, *.n) 

!■/») } 

~~ Ay )(».; + J.fc.n) — ( j+n) j ~ i.n) j = 0 

(23) 


— ^ ^ A + —g + j ) <p r * +(7 + 1 ) 0 . 
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(4> y ) and (dj) are always calculated from a centered expression. However, 
the approximation of (<p T ) strongly depends on whether locally, at a point, 
the flow is subsonic, supersonic, sonic or it is a shock point. Extending 
References 33, 34 methodologies to our case and using Eq. (21), a centered 
approximation and a backward expression are given for u": 

u U.j,k. n ) = A’ + f(j,k,n) + -g(j,k,n) 

1 7 

7 + 1 

+ -T + 1 .j.k.n) - o(i - 1 ,j.k,n)) 


u (>.j,k,n) = A' + f{j. k. n ) + ^g(j. k , u) 


+ k, n) - <p{i - 2,j, k,n)) (24) 


The local type of the flow is determined by the following table: 33,3,4 


condition 

a’ 1 -' 

u mb 

local flow is 

1 

< 0 

< 0 

subsonic 

2 

> 0 

> 0 

supersonic 

3 

> 0 

< 0 

a sonic point 

4 

< 0 

> 0 

a shock point 


Table 1 


Eq. (23) is developed in a specific form according to the local type of the 
flow. When the flow is locally subsonic, an elliptic difference form is used: 
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G(i,j, k , n) 


Ax* A + ~ o( ' 1 ~ L J'*' n ~ !)) 


When the flow 
G(i,j,k,n) = 


1 




( g(j,k,n ) — g(j, k, n - 1)) 


+ ( a ’ + fij- k. n) + -g(j. k. n) 


+ ( °( 1 + I'J' k- »)-<?(«— 1, j, k, n))| ■ 

ojj + hjjc^n) — 2o( * , j, fc, n ) + d>( z — 1 , j, k, n ) 

( Ax* ) 2 

1 

+ 1 -k, n) - 2 o(i,j,k.n) + q(i.j - 1 ,k,n)) 

~ (A-) 2 + !- n ) ~ 2o(i.j,k.n) + <p(i,],k - l.n)) = 0. 

( 25 ) 


is locally supersonic, a hyperbolic difference form is used: 


2 

Aj*"a ~ ^ 1 ~ 1 A' k. n)) - 


2 

“7“ O z -(i-j'k, n - 1) + 


1 

7 A <* 


(g{ij,k,n) - g{i,j,k,n - 1)) 
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+ { A + f{j* k. n ) H — g(j*k.n) 


{0(1. j, k, n) - o(i - 2 .j,fc.n))| ■ 

!h ’ll ~ Ifli ~ fc' n ) 4- <£(i - 2, j, fc, n) 
( Ax' ) 2 


~ ( Ay) 2 ^*' ^ ^ n ) ~ 2o{i.j. fc. n) + <p(i.) - 1 . fc. n)) 

- (A , )2 MU-* + !- n ) “ 2d(i.j,k,n) + <p{i,j,k - 1 ,n)) = 0. 

(26) 

When the flow is locally sonic, {-/-hi )<z> r . + 7 \ + / -t- = 0 . Then the sonic 

point difference form is: 

G(i,j,k,n) = ) 2 (c>(i, j + l.k.n) - 2o(i,j , k, n) + o(i,j - 1 ,k,n)) 

+ 2 (<p(hj\ A- + l.n) - 2o(i,j,k.n) + <p(i,;,fc - l,n)) 

1 / 7 — 1 

— 7 r— 7 — ( y(y. A\ n ) — y(y. A.\ n — 1 )) 

(7 + 1) A f* \ 7 

-2(f(j.k.n)- f(j,k.n- l))j =0. (27) 

When locally there is a shock point, a shock point difference operator is 
used where the flux ahead of the shock may be approximated by a backward 
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formula and the flux behind the shock bv a centered formula. 


G{i, j, k, n) = + l.j. k. n ) - 6 {i - 1. j. k. n)) 

2 1 

~ AT °*' {t -J- k - n ~ l) + ~ 9 tiik*n- 1)) 

+ |/v + f(j.k.n) + ~9(j, k,n) 

7 + 1 

+ ^-^ — 7! o( 1 + L j,k\ ti) — o(i.j.k.n) + <?( i — 1. j,k.n) 

+ o{i - 2.), k. n))| • 

g(t + I,;, L n) - 0(2,7. fr, n) - <?(i - 1. j, fc. n) + <^(f - 2, j, fc, n) 

(Ax') 2 

~ ( ^ ^ ( o( 2, 7 + l.k.n) - 2 o(i.j,k,n) + 0(2,7 - l,fc,n)) - 

A'.")i ( °( '• ./• + 1 - ” > - 2 o{i.j.k.n) + <t>{i.j,k - 1, n)) = 0. 

(28) 


Stating from initial conditions that describe a given shock wave in the 
space for t = 0 (or n = 0), and given the functions f(y, z, 0 and g(y , z, t ), 
equations (25) through (28) can be applied for n = 1 at any mesh point 
according to Table 1. They can be solved by an iterative point or line - 
or plane - relaxation algorithm until at any point max |G(t,7, k, 1 )| < 6 
where 6 is a given small tolerance of convergence. Then d x .(i, 7, fc, 1) can 
be calculated at any mesh point and the process is restarted for the next 
time step. In this way the shock motion and pressure field behind it can 
be integrated in space and time and the effect of various deterministic and 
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random fluctuations / and g can be studied. Numerical calculations of 
various examples are underway and will be presented in a future report. 


5 Conclusions 


The review of the theoretical studies of the interaction of shock waves with 
free stream vortical flows or turbulence shows that this complex nonlinear 
interaction is still an open problem to analyze. The analysis of the lin- 
earized problem of the interaction of a weak shock with relatively small 
disturbances shows that it is an invalid approach when the perturbations 
are of the order of shock strength. However, in a coordinate system mov- 
ing with the basic weak shock, the problem may fit the transonic theory 
framework. 

A new transonic small-disturbance model has been developed where a 
rescaling of the axial coordinate and time has been considered to capture the 
basic nonlinear effects that occur in the flow across the shock. This model 
results in two alternative approaches: (1) an equation for solving a velocity 
potential function that is described by an extended version of the classic 
small-disturbance equation for unsteady transonic aerodynamics; 33 and (2) 
a nonlinear equation to describe the pressure field that is similar to the 
model equation recently presented by Pierce 1, using logical considerations 
only. This equation also has extended form of the classic equation that 
describes the propagation of nonlinear sound beams with narrow angular 
spectrum. 30,31 

Both approaches show that diffraction effects, nonlinear steepening, fo- 
cusing and caustic effects and random induced turbulence fluctuations in- 
teract simultaneously to determine the development of the shock wave in 
space and time and the pressure field behind it. Turbulence fluctuations 
tend to change the local speed of sound in the flow across the shock and 
through this effect to reduce or magnify the strength of the basic shock. 

A finite difference scheme that uses Murman and Cole 3 ' 1 finite-difference 
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techniques for solving mixed- type elliptic- hyperbolic flows with shock waves 
has also been presented. Numerical calculations of the interaction of shock 
waves with various deterministic and random fluctuations will be presented 
in a future report. We will also look for analytical methods to identify the 
basic relations between given turbulence properties and the development 
of the shock waves and the pressure field behind it. It is expected to 
find solutions that will describe both peaked or rounded or messy pressure 
signatures as were recorded in experiments. We also intend to extend the 
model to include humidity and winds effects that are also known to have a 
significant effect on sonic boom pressure peaks and rise times. 
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